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Artificial Protein Cavities as Specific Llgand-binding 
Templates; Characterization of an Engineered 
Heterocyclic Cation-binding Site that Preserves the 
Evolved Specificity of the Parent Protein 

Rabi A. Musah 1 , Gerard M, Jensen 1 , Steven W. Bunte 2 , 
Robin J, Rosenfeld 1 and David B. Goodin 1 * 

Cav ity complementation has been observed in many proteins, where an 
appropriate stnail rt ol s I to a cavity-forming mutant. Here, the 
binding of compounds to the W191G cavity mutant of cytochrome c per- 
oxidase is characterized by X-ray crystallography and binding thermodyn- 
amics. Unlike cavities created by removal of hydrophobic side-chains, the 
W191 G cavity does not bind neutral or hydrophobic compounds, but dis- 
plays a strong spe< fi i ! si in with the role of 
the protein to stabilize a tryptophan radical at this site, Ligand dissociation 
constants for the profonated canonic state ranged from 6 fiM for 2-amino- 
5-methyithiazole to 1 mM for neutral ligands, and binding was associated 
with a large enthalpy-entropy compensation. X-ray structures show that 
each of 18 compounds with binding behavior bind specifically within the 
artificial cavity and not elsewhere in the protein. The compounds make 
multiple hydrogen bonds to the cavity walls using a subset of the inter- 
actions seen between the protein and solvent in the absence of bgand. For 
all ligands, every atom that is capable of making a hydrogen bond does SO 
with either protein or solvent. The most often sees 1 is to Asp235, 
and most compounds bind with a specific orientation that is defined by 
their ability to interact with this residue. Four of the ligands do not have 
conventional hydrogen bonding atoms, but were nevertheless observed to 
orient their most polar CH bond towards Asp235. Two of the larger ligands 
induce disorder in a surface loop between Prol90 and Asnl9'5 that has 
been identified as a mobile gate to cavity access. Despite the predominance 
of hydrogen bonding and electrostatic interactions, the small variation in 
observed binding free energies were not correlated readily with the 
strength, type or number of hydrogen bonds or with calculated electrostatic 
. ergii . done Thus as with naturally occurring binding sites, affinities to 
W191G are likely to be due to a subtle balance of polar, non-polar, and sol- 
vation terms. These studies demonstrate how cavity complementation and 
judicious choice of site can be used to produce a protein template with an 
unusual ligand-binding specificity. 
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Introduction 

Several recent studies have demonstrated the 
feasibility of creating artificial cavities within pro- 
tein structures that possess well-defined affinity 
and specificity for binding small molecules. 1 ** 
Further characterization of the properties of these 
artificial binding sites will enhance our under- 
standing of baste features of protein folding and 
stability. For e<ampte, this approach could provide 
an important advance over mutational analysis of 
protein-ligand interactions, as a large number of 
ligands can be used that vary in bydrophobicity, 
size or shape in a much more subtle -way than is 
possible with amino acid substitutions. Engineered 
cavities may also have practical uses, such as in 
the development of specific biosensors, where a 
designed protein template is used to specifically 
detect molecules of interest. Finally, a well-defined 
system consisting of art engineered protein receptor 
and an array of ligands thai bind with known, 
structure and energetics could aid in the develop- 
ment of toots for drug design. Extraordinary efforts 
are currently underway in the pharmaceutical 
industry to identify and optimize small molecule 
binding to specific substrate or effector binding 
sites.* Alternative approaches have been reported 
in which small molecules are identified that inter- 
rupt specific protein-protein interactions.*'' For 
both of these approaches, it is important that bind- 
ing potential is localized to smalt regions of the 
interface. However, it is not entirely clear how uni- 
versal factors such as hydrophobicity, hydrogen 
bonding, solvation, electrostatics, and protein con- 
formational dynamics combine to distinguish 
evolved natural binding sites from more common 
surface clefts. 8 Even more uncertain is five degree 
to which cavities created by mutation will resemble 
natural binding sites with respect to these proper- 
ties. 

In parallel with experimental strategies, compu- 
tational approaches are providing ever more accu- 
rate predictions of prot< n sta >i!it> Lgand-binding 
affinity and conformation, Developments in com- 
putational methods such as free -energy 
perturbation 9 * 10 and X-dyamics n and in more tra- 
ditional docking 12 and molecular mechanics 
simulations" have advanced to the point of pro- 
viding useful tools for screening drug candidates. 
These methods have been primarily developed and 
tested using data for li ig to naturally 

evolved sites ini 'm i enzymes WhUe such 
sites have provided a particularly appropriate 
benchmark for computational method develop- 
ment, tests of their perfonnance with ligands 
bound to artificial protein cavities may be of sig- 



nificant interes t , i sites may display different 
compositions - ■ wc of polarity, 

conformational flexibility and solvation character- 
istics. Irs addition, the lower affinity and specificity 
expected for ligand binding to artificial cavities 
should provide unique and challenging tests of 
o mputatiooal pn ii ri< ns 

Detailed characterization of the interactions 
between ligand and protein is necessary before 
ligand bfoding to artificial cavities may be usefully 
compared to their natural counterparts. In previous 
studies of artificial prot< i ■ - tw o rather dis- 
tinct types have been described. Largely hydro- 
phobic huriccl cavities created in phage ?4 
lysozyme 5 - 114 ' 15 generally do not appear to be 
speuhtalh 1 a tspt r> a mixture of rigid 

and eonformationally mobile characteristics. Such 
cavities appear to bind hydrophobic ligands with a 
specificity that is determined by a combination of 
packing, desolvation, and conformational entropy 
factors. On the other hand, several polar buried 
cavities have been created in the enzyme cyto- 
chrome c peroxidase (CCP), each of which speetfi- 
i i " ! One such cavity, 

W191G, created by deletion of foe side-chain of a 
tryptophan radical cation center is shown in 
Figure 1. In contrast to the hydrophobic cavities of 
14 lysozyme, several molecules of solvent and a 
monovalent cation occupy the W191G cavity 
mutant of CCP in the absence of ligand, and the 
cavity was observed to bind small cationic hetero- 
cyclic compounds. 3 '* in many respects, the inter- 
actions seen in the polar W191G CCP cavity are 
more reminiscent of natural substrate-binding sites, 
with the potential for specific hydrogen bonding 
and electrostatic interactions. It has not been dear, 
however, how promiscuous the W191G cavity is 
with respect to ligand variation, whether more 
than one mode of binding is utilized, what range 
of affinities can be achieved, or how much the sur- 
rounding protein structure responds to ligand vari- 
ation. Here, we report the structural and 
thermodyrsamic characterization of a variety of 
ligands bound to this cavity to address these ques- 
tions. 



Results 

Specificity and thermodynamics of 
ligand binding 

Initially, a wide variety of small organic mol- 
ecules wese screened ! tity to bind to 
W191G. Previous studies have shown thai displa- 
cement of the solvent within the W191G cavity by 
ligands such as imidazole results in a subtle but 
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characteristic perturbation of the heme Soret 
absorption band. 3 Thus, samples of W191G were 
screened by titration against potential ligands 
selected primarily by their ability to fit into {he 
known volume of the W191G cavity. As controls, 
identical titrations were carried out against WT 
CCP and the W191F mutant, the latter of which 
has heme Soret absorption features similar to those 
of VV191G and can be considered to be a model for 
the cavity that is irreversibly occupied by the Phe 
side-chain.^* 13 These features of the control 
were particularly useful in discriminating binding 
from non-binding ligands. Those compounds that 
produced the characteristic optical perturbation, in 
W191G but not W191F were designated as poten- 
tially binding. As a result of this analysis it became 
apparent that previous conclusions^ about the 
determinants for ligands to W191G were accurate; 
namely, that ti - pecific I •- n\i!l hetero- 
cyclic cations. A number of these, such as imida- 
zoles, exist in equilibrium between neutral and 
pmtonated eationic forms with physiologically 
accessible pK, values. In each such case tested, 
ligands appeared to bind more tightly to W191G 
below their respective pK, value than above U 
(data not shown). Based on these results, a second 
round of candi > tocusmg on 

shape-selected compounds with physiologically 
accessible eationic forms to produce a list of 
approximately 20 compounds that appear to bind 
specifically to W191G. From these screens, the fol- 
lowing: trends were readily apparent: (i) all tested 



heterocyclic candidates that are small enough to fit 
into the cavity and are eationic under the exper- 
imental conditions elicited the optical binding 
response; {») neutral compounds, such as indole, 
or compounds with a pK a well below the pH of 
the titration (typically pi t 4.5), such as cyanopyri- 
dine or tetrazde, did not bind; and (hi) eationic 
compounds that are non-planar or that would 
require protein conformational change to alter the 
cavity dimensions, such as rrimethylarnine and 
pyrroline, gave no evidence in support of binding. 

Dissociation constants were measured for com- 
pounds identified to bind W191G, using two inde- 
pendent methods. Measurement of the 
perturbation in the heme Soret absorption was 
used to obtain dissociation constants in titrations 
of W191G with ligands, and for selected com- 
pounds, binding param 1 > \ i btained by bo 
thermal titration calorirnetry (ITC). These data are 
presented in Table I, where entries mat include 
enthalpies and entropies were derived from ITC 
measurements and the remaining entries were 
from optical titrations. In agreement with the initial 
screens and with previously published data, 3AM - w 
it was clear that only compounds that were predo- 
minantly eationic at the measurement pH of 4.5 
bound to the protein, while those that were neutral 
under the experimental conditions (i.e. pK, below 
4.5), were observed not to bind. Thus, ail data pre- 
sented in Table 1 have been corrected for the pA' a 
of the compound by taking into account the con- 
centration m solution of the eationic form, so that 



Ligand Binding to an Artificial Protein Cavity 



Table i < \ r j i < - ! g u 1 1 i n, „■ 



2a5mt 2-Am.uKvf!-ffi«hyUhi;i^o1f 

l2dmi 
2ai 

4ap 4-Amirir>pyndw 

3Ap 3-AminopyrkUne 

J mim "t -Methy timki ii xole 
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these values represent intrinsic binding affinities of 
the canonic forms. Overall, the agreement between 
the optica'! and calorimetrtc data is good. Dis- 
sociation constats range from 6 jtM for 2-amino-5- 
methylthia<s.oit < appro* twtdy 1 raM in 
the case of 2,3,4-trimethylrbjazale (234tmt) for 
those ligands shown to bind. This represents a 
range of 3 kcal/mol (1 cal = 4.184 j) in binding 
fn.t< energy. Both the binding enthalpies and entro- 
pies were observed to he negative, as is -usually 
observed in ligand-pratein binding/" One signifi- 
cant observation is that there is a much greater 



variation in binding enthalpy than for the total free 
energy, which implies a significant entropy-enfhat- 
py compensation, as shown in Figure 2. 



Structural characterization of ligand binding 

Ligand binding to the W19IG cavity was charac- 
terized by X-ray diffraction using crystals soaked 
in mother iiqs . tain > i is Uganda. Pro- 
cessed data sets (see Table 2) for ligand-soaked 
crystals were matched 1 < , \y det^nmtud 

structures of tinsoaked WI91G crystals of the 
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Figure 2. A plot of the thcnnodyrwrnic parameters 
determined for ligand binding to the W191G cavity. 
Estimated erroi bat town in the bottom graph for one 

- '■- btained horn four repealed 

measurements. This plot illustrates that AH values vary 
by more titan the total free energy, AG, Thus, a ssrge 
entropy-enthaipy compens Boft i bserved as shewn in 
the above graph. 



appropriate crystal form,' 1 Ligand binding tvpicallv 
displaces the 'four well-ordered water molecules 
and potassium ion observed in the cavity in the 
absence of ligand, as shown in Figure 3(a). Thus, 
direct F s .,, v , - F ummV difference maps were often 
not easily interprelable because both positive and 
negative density features are superimposed. 
Instead, a series of — F e omit maps were cre- 
ated, where f^.,. are the observed structure factors 
of the ligand-soaked crystal and F,_ are the struc- 
ture factors derived from a model for unsoaked 
W191G that was refined without any atoms or 
water molecules included in the cavity. Such omit 
maps represent the electron density of the contents 
of the cavity under the soaking conditions. These 
omit maps, contoured at 4c and Set, are shown in 
Figure 30>)-(s) f< \e mds bat were desig- 

nated as potential f binders.* Seariy interpret- 
able density is observed in each case, indicating 
occupation by the sp< > Sij i fi i fly, 

no difference density above about 2c was observed 
with any ligand in any other region of the struc- 
ture remote from the W191G cavity, clearly 
demonstrating that alt binding interactions are 
is in spite of 

the fact that is might be capable 

of coordinating the heme iron by occupying the 
naturally occurring channel that' is distal to the 



heme (partially visible in Figure 1). In addition, 
test soaks with ligands that were designated as 
non-binders in the initial screens, for example tetra- 
zote and 3~anhnopyra<tole, showed only the pat- 
tern associated . it • t in the unsoaked cavity 
(Figure 3(a)), 

fa almost all cases, it was possible to assign the 
orientation of h^am si in the on t> Fur thnv 
ligands, including 1,2-dimcthylimidazoie {i2dmi), 
aniline, and 2-etl I la ok (2eim), the shape of 
the electron density envelope was sufficient to 
assign ligand orientation. For others containing sul- 
fur atoms in the ring, the position of the sulfur 
atom was evident in the electron density at a high- 
er contour level, as shown in red in Figure 3. From, 
this additional information, it was possible to 
assign ligand orientation for 2-aminotfuazole (2at) ; 
2-3mino4-methylthiazole (2a4mt), 3- methyl fhia-- 
7,ole (3mt) XI U t i 1th 1 . 04Stmri and 
23,4-trim.ethyfmiazole (234tmt). Finally by combin- 
ing the above information with inferences derived 
from hydrogen bonding interactions, as discussed 
beiow, orientations of 2-am irto~5- methy ! th ia zole 
(2a5mt), 3-aminopyridine (3ap), 1-methylimida^ole 
(Imim), 2-aminopyridine (2ap), 3,4-dimethy- 
IthiazoSe (34dmt), 2~methyhmidazo)e (2miivi), 
imidazotU-ajpyridfoe (imp),'indotine, and 1-viny- 
Itmidazole (Ivjm) can be assigned unambiguously, 
Of the IS ligands, 4-aminopyridine (4ap) cannot be 
assigned to a specific orientation, and lvim 
,'1 >K >d density to assign 
the conformation of the vinyl group. These 
observations may result from lack of resolution, 
structural landmarks, or true conformational het- 
erogeneity within the cavity. Nevertheless, reason- 
able estimates of the most likely orientation for 
these compounds were possible. Based on these 
assignments, models for the ligands were manually 
placed into the electron density as shown in 
Fig ,ie 3(faHs) 

With few exceptions ligand binding to W141G 
displaces the solvent of the ligand-f ree cavity with- 
out changes in the cavity wall Previous studies 
have shown thai W191G in the absence of ligand 
contains four water molecules, W401 -W404, and 
one monovalent cation, either K* or Na + that are 
well ordered in the cavity and form a network of 
hydrogen bonds between themselves and with 
polar groups on the cavity walls. *<* These inter- 
actions are listed trt the first row of Table 3. One 
wafer molecule, W30S shown in Figure 3a, was nest 
introduced by the cavity but is observed in wild- 
type (Wf) CCP and is "retained in W1MG. With 
me exception of 2a t, 3ap, and aniline, each of the 
ligands displaces all of the solvent atoms intro- 
duced by cavity formation, without displacing the 
conserved W308, For 2at, 3ap, and aniline, one of 
the cavity water moki ' . emams in the 
cavity in addition to the ligand. While most 
ligands bind to VV191G without inducing changes 
in the dimensions or shape of the cavity walls, two 
of the ligands, ktdoline and imp, appear to have 
induced an alternate protein conformation in the 




Figure 3 (legend shoum on page 8S2) 




Figiwe 3 {kgend shown on page 852) 
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cavity. For these two structures, significantly 
reduced elect . . , seen in the protein for 
residues 190495, indicating that this region of the 
structure has become mobile or is occupying mul- 
tiple conformations. 

Several modes of hydrogen bonding interactions 
are observed in the crystal structures of Figure 'A. 
Hydrogen bonds inferred from the assigned tigand 
orientations are listed in Table 3, where the rows 
are sorted in order of the binding affinity, and the 
columns ordered by the number of total inter- 
actions. The most common hydrogen bond is to 
the earboxylate side-chain of Asp235. This inter- 
action occurs for 12 of the 18 ligands and for two 
of them, 2a5mt and 2ap A.sj 13! jppears to make 
a bifurcated hydrogen bond to two atoms in the 
ligand. The second most often observed interaction 
between ligand and cavity, which occurs in six of 
the ligands, is a hydrogen bond to the conserved 
water molecule of WT CCP, W308. In addition, a 
few of the ligands make hydrogen bonds to the 
backbone carbonyl atom of L177, H175, and M23G. 
Each of the three ligands, 2at ( 3ap, and aniline, 
wl-uch retain cavity water W401 also form hydro- 
gen bonds with it. Pour of the ligands, 34dmt, 
343tmt, 3mi, and 234tmt, are methylthiaozies that 
are cations without traditional hydrogen bonding 
capability. However, these four compounds orient 
their most acidic CH proton (that adjacent to the 
endocycMc ring nitrogen) so that it is directed 
towards the \ ; >p23 x>xy late group. Finally, it 
is noted that every atom in each of the 18 ligands 
that is capable of forming a hydrogen bond does 
so with some atom in the protein or included sot- 
vent. 



Discussion 

The foremost conclusions from this study con- 
cern the unique specificity for ligand binding to 
the artificial cavity of v\ 19 li . Several cavities have 
been introduced into proteins by deletion of amino 
acid side-chains, and it is often observed that these 
cavities bind small n ol at m nic the prop- 

erties of the removed structure. 1 ""'" 32 Titus, artifi- 
cial cavities can be utilte-ed as templates for 
binding behavior. For example, several cases have 
been reported in which amines rescue function of 
lysine residue* that have been replaced by smaller 
side-chains. 15 Other artificial cavities that have 
been created by removal of side-chains in the 
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interior of a pr • t ne rally been non-polar 

in nature, reflecting the properties of the buried 
residues removed to create them."' 14 ' 5 ' These 
non-polar cavities do not appear to be specifically 
salvated, and bind h\ I b gauds without 
making specific in i 1 ' moderation of 
hydrogen bonding and electrostatic effects are 
minimized and the crucial concerns are with 
hydrophobic interactions, packing and confor- 
mahonal entropy. Uo>m i sty introduced 
by the W191G mutation of CCP is quite different 
in these respects. The Trpl91 side-chain, while 
itself quite hydrophobic, is a special site in the 
enzyme that is oxidized reversibiy to a stable rad- 
kal'catiort form as a part of the enzymatic mechan- 
ism.*'* Not surprisingly, the structure surrounding 
this sste has been evolved specifically to stabilize 
the cation form. 4,27 "' 30 It is remarkable that the 
ligand-binding specificity of the W191G cavity 
faithfully reflects this evolved property, and that 
all of the 20 compounds identified to date that 
bind to W191G are cationic heteroeycies that bind 
only to the cavity by making an array of hydrogen 
bonding interactions. 

Subsets of the same interactions that are believed 
to stabilize the naturally occurring Trpl91 cation of 
WT CCP aw utilized to define the specificity deter- 
minant of VY191G. We conclude that electrostatic 
interactions with Asp235 and several backbone car- 
bonyl atoms define the cation requirement/' 1 * 5 " 32 
while specific hydrogen bonding interactions to 
these groups determine the ligand orientation 
within the cavity. An important observation is that 
the highly ordered solvent observed in the W191G 
cavity in the absence of ligand forms an ideal tem- 
plate for ligands, because all potential hydrogen 
bonds to the cavity walls are satisfied. While none 
of the discovered ligands makes all of these poten- 
tial hydrogen bonds simultaneously, no new 
hydrogen bonds an uti ' i 1 ; ij ands that neie 
not also seen with solvent alone. Of the six differ- 
ent hydrogen bonds made to ligands, Asp235 
clearly dominates and appears to determine ligand 
orientation- This is seen from comparisons of the 
isosteric ligands 2ap, 3ap, 4ap, and aniline, where 
each ligand adopts a different orientation in order 
to permit a hydrogen bond with Asp235. Similarly, 
2eim and Ivim orient their substiruents differently 
within the cavity in order to preserve interaction 
with Asp235. Finally, 34dmt, 345fanl, 3mt, and 
234tmt, are cationic methylthiaozles without tra- 



Figwte 3. A series .if stereo electron density omit maps showing the contents of the W191G cavity after soaking 

i-waS u a u f th<: cavity tor a crystal without swaking in 
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TaMe S, Hs< >^hh Interact r f een file W191G «avity and bound ligands 
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ditional hydrogen bonding capability. However, 
each of these four compounds orients its most 
acidic CH proton {that adjacent to the endocyclic 
ring nitrogen) so that it is directed towards the 
Asp235 carboxylate group. This interaction, dis- 
cussed previously in terms of an unconventional 
CH to O hydrogen bond," illustrates the remark- 
able ability of the ligand-cavity interactions to 
mimic hydrogen bonding patterns specified by the 
template of the removed tryptophan side-chain, 

Calorimetry data show that ligand binding to 
W191G is enthalpy-driven, but the overall binding 
free energy is relatively weak, does not vary 
greatly, and exhibits large enthalpy-entropy com- 
pensation. These factors make empirical corre- 
lations of the relative binding energies difficult to 
predict without detailed calculations. Entropy- 
enthalpy compensation is commonly observed in 
aqueous solution, and is usually characterized by 
slopes near unity, 93 as observed in this case. This 
effect is often attributed to a trade-off m which a 
stronger enthapHc interaction results in a compen- 
sating loss of entropy due to motional restriction. H 
It is expected that the solvation energies of the 
ligands used in this study will be quite similar 
Indeed, calculations using MEAO^' (data not 
shown) indicate a variation of only about 20% in 
the solvation energies of these compounds. Simi- 
larly, no readily apparent relationship exists 
between the number or type of hydrogen bonds 
made to the Iigand and the relative thermodyn- 
amic parameters. This is not unexpected, as these 
very polar compounds will be fully hydrogen 
bonded in solution, and have also found fully sat- 
isfied hydrogen bond , hit the cavity. 
Thus, while hydrogen bonding interactions are 
je absolute interaction 
enthalpies with the cavity and may dictate ligand 
orientation, differences between bound and 



unbound interactions will be similar for each com- 
pound. A more sophisticated attempt was made to 
con-elate electrostatic interactions with binding 
enthalpies. As described previously, electrostatic 
calculations using POLARIS and the FOLD {pro- 
tein dipoles l angevin dipoles) method 3 *" 3 * were 
performed to estimate the electrostatic interactions 
between the ligand and the surrounding protein. 
These calculations {not shown) also failed to 
■ cou it foi va 5i? ions in tht 1 rved enthalpies or 
free energies. The small difference (—3 kcal/mol) 
in binding free energy between the tightest and 
weakest of these ligands most likely reflects a com- 
bination of subtle interactions including electro- 
statics, hydrogen bonding, desolvatkm entropies, 
and packing interactions that must each be 
accounted for very accurately in order to predict 
relative affinities. Previously successful appli- 
cations of the FOLD method in comparing the rela- 
tive binding of 234tmt and 345tmt 5 * or in the 
comparison of the relative stabilization of the 
tryptophan cation radical in CCP relative to ascor- 
bate peroxidase 32 likely resulted from the cancella- 
tion of elements such as solvation and non-polar 
interactions that are not well modeled by the 
specific implementation of the PDLD method used 
in these studies. 

As with the previously studied cavities in T4 
lysozyme mutants, the W191G cavity contains 
regions that appear structurally rigid, and others 
that are subject to movement. It might be expected 
that artificially cj ted ities will undergo some 
collapse, or be more confotmationally mobile than 
natural binding sites. However, a detailed analysis 
of cavity collapse has been reported for T4 lyso- 
zyme mutants, 15 and it is remarkable that the cav- 
ities that have been characterized to date generally 
appear quite rigid. One such cavity, the L99A cav- 
ity of 14 lysozyme, was observed to have a region 
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that is structurally rigid and another region that 
deforms in response to - scking T^ractiora with 
ligands. M In the case of the WK'IC cavity of CCP, 
very little, if any, collapse or movement in the 
dimen-kjns i t seen upon binding for 

most of the Jigands. However, the cavity is inac- 
cessible to solvent in the average crystal structure, 
and thu*> the relatively :ap>d excha 
requires a significant role for conformational 
dynamics in order to i 1 > ess ir Indeed, 

when bemunidazole is bound to W191G, a large 
bop rearrangement between Ptt)19f) and Asnl95 
has been observed to give an open channel confor- 
mation, providing a tantalizing view of the path- 
way for ligand access.. 17 It is of interest that the 
two largest" Ugands of this study, imp and indoline, 
show significaniiy reduced electron density for the 
protein at thfc rt r i uggestii\g confor- 
mational heterogeneity of this loop. Thus, this cav- 
ity appears to contain a specific hinged gate with 
open and closed conformations to allow ligand 
access, while the remaining walls of the cavity 
appear quite rigid and resistant to either collapse 
or induced fit binding behavior, 

The specif i ting properties of the 

W191G cavity may help establish practical uses for 
engineered cavities, for example, a number of fac- 
tors make the W191G cavity an ideal test bed for 
use in the evaluation and development of compu- 
tational tools for drug design. The rigidity of the 
cavity wall upon bin 1 i nphiies the 

task of accounting for cavity collapse or adjustment 
in the protein structure. The polar nature of the 
cavity and its Hgands provide a template that is 
similar to naturally occurring enzyme sites. Tire 
small variation in binding free energies provides a 
challenging test for calculations of the delicately 
balanced forces involved. Finally, the crystal struc- 
tures of a complete series of Ugands bound to the 
same site along with the thermodynamic par- 
ameters provides a well-defined framework for 
such calculations. Computational studies of ligand- 
binding energetics have already made use of the 
W191.G cavity mutant for method development. 
The.se studies have resulted in improved compu- 
tational methodologies and have provided specific- 
predictions concerning potential conformational 
heterogeneities for certain Hgands that are testable 
by further experiments/ 57 ' 38 



Experimental methods 

Protein expression and purification 

W191.G ap< instituted with 

heme, and purifi led enzyme was 

recryMalliaed twice from distilled water and stored as a 
crystal suspension at 77 K. Before use, crystals were 
iissolved in the 

appropriate buffer. W19'iC protein concentrations were 
M~ cm - ' as determined 
from the pyridii < c nogen assa> 



Ugand-binditig measurements 

Two methods were used t I- terUe ligand bind- 
ing to the WW i t > t> b ition constants 
were obtained by optically detected titrations as 
described 1 If est ip' n a <wall 

perturbation of the hern ■ >rbance that t$ 
observed when ligands displace the solvent in the 
W191G cavity lb mea rei nts ligand stock sol- 

utions were prep; e I mol, except for 

indoline, imida i incline which 

were prepared in 95 % ethanol. Prior to use, each stock 
solution was adjusted to the same pH as the protein buf- 
fer with HjPCV Dissociation constants were determined 
from Scatchard t . the difference absorbance 

of she Sorel maximum, assuming one binding site per 
protein molecuk For selected ligai s more extensive 
characterization oi the Hgai 1-blnd g therrood; mi 
was obtained by isothermal titration caiorbtnetry (TIC). 
These meisurements. wi c . < M< I titra ion 

akititiwte from Microcal, Inc. at 25 C in lOfi mM Bis- 
Tris propane at pH 4, c . - of known con- 

centration were titrated with known concentrations of 
hgand (typh il t t on relative to 

p.t tea .ii led 5 i'in|ecti >r«s from a hX> uf spin- 
ning syringe (400 rpm) at intervals of four minutes, 
Titrations were extended past the end-point to allow 
subtraction of the heat of ligand dilution from the injec- 
tion peaks Data itaij is ■> 3 performed using ORIGIN 
software customized for FTC analysis by fitting the data 
to a single-site binding isotherm to obtain values for the 
nd assodatfo ' r 1 thalf nd t!i 
number of ligand sites. Binding free energies were caicu- 

t f x ' i -state and bind- 

ing entropies were calculated as (AG - &H)/T. With the 
- V j i of t-dimetby StfuaEoie and 2,3,4- and 3,4,S-trh 
tnelhylthiazoles, which were synthesized by known 
1 I t i i i h. As 

the results sh > y m a cattonic 

form (see Results), each of the reported values for bind- 
i >s been torrccted v, here 
appropriate for the pK, of the ligand to reflect the intrin- 
sic affinity of the cationic form. 

X-ray crystallograprtic data collection and analysis 

Crystals of W)5 t wn fo < tw dai 1 ) culler 

don by vapor diffusion at "W'C from sitting drops of 
CCP in 8.5% fv/v'J 2-methvi-2 ) 4-pentanediot (MPD), 
200 mM Mes , F H $.0), against *er yk of 25% MPD, 
as described. 5 C < tab i - a ed a' artificial mother 
liquor containing 30-50 mM ligand for 60 mirvutes before 
mounting md d . i < Hffracti n data were col 
lected at 15 *< , K li I , t the rotating 
anode of a Siemens SKA X-ray generator using a 
Siemens area detector. Data were indexed and integrated 
using the XF.N 1 1 yens analyzed 

by difference Pouri i teehnkjui - > sing the Scripps Xtal- 
View scrftware.*' Models fot the iigsmds were obtained 
by geometry optimisation using density tunctionai 
methtxfs implemented jn C n-9<J is desenbed w 

Protein Data Bank accession cod«s 

The crystallograpluc coordinates for tire structures 
presented in this work have been deposited with the 
RCSfl ViiAem Data flank (http://www nrslmrg) with 
accession codes laeq, laej, laee, laeo, laef, laeg, laeb, 
iaed, laen, laeh, laek, and laem. 
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